Constraints on the three-fluid model of curvaton decay 



T. MultamakiQ J. SainioQ and I. ViljfH 

Department of Physics, University of Turku, FIN-20014 Turku, FINLAND 
(Dated: February 2, 2008) 

A three fluid system describing the decay of the curvaton is studied by numerical and analytical 
means. We place constraints on the allowed interaction strengths between the fluids and initial cur- 
vaton density by requiring that the curvaton decays before nucleosynthesis while nucleosynthesis, 
radiation-matter equality and decoupling occur at correct temperatures. We find that with a con- 
tinuous, time-independent interaction, a small initial curvaton density is naturally preferred along 
with a low reheating temperature. Allowing for a time-dependent interaction, this constraint can be 
relaxed. In both cases, a purely adiabatic final state can be generated, but not without fine-tuning. 
Unlike in the two fluid system, the time-dependent interactions are found to have a small effect 
on the curvature perturbation itself due to the different nature of the system. The presence of 
non-gaussianity in the model is discussed. 
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The problem of determining the evolution of large scale perturbations in a background of a multi-component fluid 
system is central in modern day cosmology [l|, 0, HJ • In such a system, the interactions between the different fluids 
, are important in determining the evolution of the curvature perturbation [i[ . Examples of interactin g flu id systems 
demonstrate the importance of such s yst ems, as most notably reheating at the end of inflation [13, ITTI Il3l H3 . [m H?} 
and the curvaton scenario fll. Il9l. I2H [2rl \2l\ . 

In contrast to any single fluid system, in a multi-component system the total curvature perturbation, £, generally 
evolves whenever the non-adiabatic pressure is non-zero, i.e. when interactions between the fluids exist. Evolution of 
I' the primordial large scale curvature perturbation can relax the underlying assumptions on the inflationary scenario. 

Therefore analysis of multi-component fluid systems may affect our view on the physical settings. In addition to the 
. curvature scenario considered in this paper, natural frameworks for such mechanism exist e.g. within a traditional 
multiple inflationary scenario [33[ or a string landscape picture [34l l35j |. Whether a given scenario can effectively 
C3 . modify the primordial spectrum depends on the exact nature of the system. 

Recent cosmic microwave surveys have also brought attention into the concept of non-gaussianity i.e. how much the 
spectrum deviates from gaussian distribution. This is especially important in multifield models of inflation including 
■ the curvaton scenario 

£N| | The mechanism how energy is transferred between the fluids can be describe d by different methods, e.g. by a 
constant interaction Q or by utilizing the so-called sudden decay approximation [l9[ H3] • In a recent paper [25[ we 
considered relaxing the assumptions behind these approximations by allowing for time dependent interactions while 
evolving the full large-scale perturbation equations. Such an approach can better model the micro physics behind a 
particular physical framework by allowing one to choose the strength and the time at which the interaction is turned 
on. In contrast, if the interaction between is modeled with a constant interaction term, the fluid begins to decay (or 
interact) when its decay width is of the order of the Hubble rate, V ~ H. Physical scenarios relevant to having time 
(and space) dependent interactions include e.g. phase transitions, multiple inflation scenarios [33| and scenarios where 
locally different decay rates of the inflaton are generated by spatially varying reheating temperature and couplings 

0QBG3. 

In the present paper we consider the curvaton scenario with time dependent interactions between curvaton and 
other fluids. During inflation the curvaton is a light scalar field that does not contribute to the expansion of the 
universe; after the inflaton field has decayed into relativistic particles the curvaton begins to oscillate and to decay 
into radiation and matter. The focus of this article is on this situation: we study how a time-dependent interaction 
affects the evolution of the curvature perturbation. We study the physically allowed parameter space by utilizing 
information from known cosmological epochs. Moreover, we calculate the amount isocurvature in terms of curvaton 
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decay widths. Finally, we also discuss how to describe non-gaussianity in the three fluid model and present the 
non-linearity factor f^L [!]• 

This paper is organized as follows. In sections II and III we present the governing equations of the background and 
perturbations in the Newtonian gauge, including discussion of non-gaussianity as described by the Jnl parameter. In 
section IV we present our results while the discussion and conclusions are presented in the following section V. 



II. PERTURBATION EQUATIONS 



We begin by presenting the equations of motion of the background variables and different density perturbations, 
where we adapt the notations and conventions of [l|, 0] as well as we consider linear scalar perturbations about a 
spatially flat Friedmann-Robertson- Walker -background in a Newtonian gauge: 

- ds 2 = g^dx^dx" = (1 + 24>)dt 2 - a(t) 2 (l - 2ijj)S ij dx i dx j , (1) 

where a is the scale factor. We have used units where 8irG/3 = 1 and c = 1. 

The background evolution is determined by the Einstein's equations, — SirGT^, and the continuity equations 
of individual fluids: 

P(a) = -BH(1 + UJ(a))P(a) + Q(a), (2) 

where lj^ — P( a )l P(a) is the equation of state and Q( a ) describes the energy transfer between different fluids. 

Perturbing the covariant continuity equations, one finds the evolution equations of the energy and pressure density 
perturbations <5/C( ) and SP( a ) on large scales |4|: 

S P(a) + 3-ff(<5p(a) + SP(a)) - 3(P(a) + P(a))^P = Q(a)<t> + 8Q(a)- (3) 

In addition, we have the G® component of perturbed Einstein equations 

ZH(i) + H<j)) = -AitGSp, (4) 

where Sp = J2 a ^P{a)- For perfect fluids, (f> = ip, and hence given the equations of state, and the interactions 

between the fluids, Q( a ), one can evolve the individual fluid perturbations along with the metric perturbation <p. 

Since gauge dependence is an important issue in perturbation analysis we have used gauge invariant perturbations 
£( a ), which equal density perturbations <5/9( a ) in uniform curvature gauge ip — 0, i.e. 

6(a) = 5p(a) + P(„)V> = -P(a)C(a), 
a 

where Q a \ is the curvature perturbation in uniform density gauge, i.e. 5p( a ) = and comma represents derivative 
with respect to the number of e-folds N) In a. 

The corresponding equations of motion for the gauge invariant quantities are thus 

C(o) = + W (a)K(a) ~ MP lnt(a) + — — + — — — 

- -3(1 + W(a) )e (a) - 35P mt{a) + - ^ ^— e (c) + - — — , 

where SP int ^ = SP( a ) ~ P'( a )^P(,a)/p[ a ) are * ne internal pressure perturbations. We have also included the possibility 
of explicit time dependence of the interaction terms. The equation of motion of the curvature perturbations can be 
derived from eqs. © and The result is 

, SSPjntja) 1 H'Q {a) 

where 6Q mt{a) = 5Q (a) - Q[ a )Sp( a )l p[ a y 
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III. EVOLUTION EQUATIONS 



The model in hand is a three fluid system with curvaton, radiation and non-relativistic matter fields denoted by 
subscriptser, 7 and m, respectively. Since the curvaton field undergoes coherent oscillations, the curvaton fluid can 
be safely estimated [24j to behave like non-relativistic matter. Thus we have ui„ = 0, u> 7 = 1/3 and u m = and 
the background equations can be written in terms of fractional densities fl a = pi a )/p for which the equations of 
motion are @: 



o 7 = fi 7 (0, 



Qa_ 

Hp 1 

-1) + — , 

Hp 

Q in 



(8) 



From the definition of £l a it can be easily seen that fl a + f2 7 + f2 m = 1, which means that one of the equations of 
motion is redundant. The interaction terms read as 

Qa = -T 1 j 1 (N)p a - T m f m {N)p„, 

q 7 = r 7 / 7 (JvK, (9) 
Qm = r m / m (iV)p (T , 

where T/ a ) denotes the strength of the interaction and functions /( a )(iV) include the explicit iV-fold (time) dependence. 

From eq. ([6]) we can finally derive the equations governing the evolution of the perturbations of curvaton-radiation- 
matter system. The equations for gauge invariant perturbations are 

c , QC r 7 / 7 (jv) + r m / m (i\o ^ r 7 /;(jv) + r TO /; t (jv) Q a $ 

C = ~3Ca Jj Ca Jj P«<t> 

?7 - 4 ?7 + jj + # P** H 2p 

U- <*U + ^ 4a + ^ /V0 -^-^ 

</>' = {— -1)0-— • 

V K H ' 2p 

Correspondingly, the equations for curvature perturbations (a in the uniform density gauge read as 

_ r„X(jVK(y, + Ca) ff' (r 7 + r m / m (jv)) ^ _ 
Cct ~ fftf, + ff ^ (C Cff) 

>/ _ r 7Pa(Ca - Cy) T -(Pa /> > v ri . s 

^~ fl^ £T p 7 U W 



(10) 



r„ 



Hp' n 



C = -EF3- fm(N)pUtm - £0 + f m (N)P«W + Cm) " _ ( C - Cff ) 



ff' Y m f m {N) Pa 



H P' m 



The evolution of the system described by eqs © and (JTU]) is then ready to be solved numerically once the initial 
conditions have been set. 

Because some of the initial conditions lead to non-physical solutions {e.g. too high reheating temperature or the 
universe might become matter dominated during nucleosynthesis) , the numerical analysis has to be performed carefully. 
In order to eliminate clearly unphysical scenarios we use physical knowledge from notable epochs of cosmology, namely 
Big Bang Nucleosynthesis (BBN), radiation-matter equality and decoupling. By fixing the temperatures at those 
times we are able to determine constraints for the energy distribution between different components of the system. 
The physical temperature scales have been set to coincide with values given in [321 ] : nucleosynthesis came about at 
temperature at least 0.1 MeV, the matter-radiation equality was reached when T = 1.0 eV and decoupling occurred 
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when T — 0.1 eV. Because p m , p a oc a 3 and p 7 oc a 4 we can derive a lower limit for the abundances of curvaton, 
CDM and radiation during nucleosynthesis. Straightforward calculation results in the limit 



n 7 
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> 10 5 . (12) 



This limit alone does not, however, guarantee that the curvaton has decayed before the nucleosynthesis. Because even 
a small contribution of curvaton during nucleosynthesis may eventually begin to dominate the system afterwards and 
therefore lead to a system with undesirable physical properties, we require that it decays rapidly enough to become 
subdominant. In practice we compel the beginning of nucleosynthesis to happen after the curvaton starts to effectively 
decay, i.e. 

A^max < N nnc , (13) 

where iVn^max is the time of maximum curvaton proportion. Setting the temperature at radiation-matter equality 
allows us to estimate the initial temperature of the system i.e. reheating temperature by T^h ~ (p~/,rh I P^,eq) x ^T eq . 
More precisely, the reheating temperatures we quote here are lower bounds since the decay of the curvaton adds to 
the radiation energy density of the system, but numerical work shows that the effect is small and in practice the initial 
temperature can well be approximated by the above formula. 



A. Non-gaussianity 

Single field inflation models usually predict the amount of non-Gaussianity to be too low to be detectable by 
future CMB-surveys. In contrast multiple scalar fields (e.g. the curvaton scenario) can lead to a clearly observable 
non-Gaussianity being therefore testable in near future. 

The non-Gaussianity parameter /jv l @, H3| is defined via gauge invariant Bardeen potential 

^ = ^ g + fNL<f 2 g , (14) 

where $ s is the gaussian part of <E>. The Bardeen potential can be expressed at the time of decoupling when the 
universe is matter dominated as $ = 3£/5. The non-Gaussianity parameter /tvl is common to express utilizing an 
additional parameter defined by 



C<7,0 



(15) 

dec 



which tells how effectively the initial curvature perturbation transfers to the total curvature perturbation. Thereby, 
using the definition of the curvature perturbation £, r-parameter, the equation of motion of p a and a second order 
estimate 8p a j ' p a — 25a /a + (Sa/a) 2 [27[, the Bardeen potential can be cast in the form 

51 + Jh7, a 4 r \ 51+ 37f^ a '°/ 

and hence /jvl — 5/(4r). Since during decoupling £ ~ Cm|Dcc and in the three-fluid curvaton model Cm|Dcc = Cx,o 
@, our value for the transfer parameter is 1 and Jnl = 5/4. As pointed out in 27] the previous estimate of Jnl is 
valid only when Jnl ~S> ^ since we are using first-order perturbation theory and we are assuming second-order terms 
$( 2 ) to be at most of order <& 2 g . Thus the three-fluid model is unable to give any limitations on the non-Gaussianity 
parameter when the linear theory is applied. 

This result differs from the non-gaussianity given in Q mainly because our /jvl is defined at the time of decoupling 
whereas in [|| the time of nucleosynthesis is used. Our analysis, however, follows the general formalism presented in 
(3l| which allows our fjyL to be compared to the observable first order Sachs- Wolfe effect [13] 

which is evaluated at the last scattering surface and where the lapse function <f> — in the Newtonian gauge. 
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IV. NUMERICAL RESULTS 



We have studied the evolution of equations (|SJ) and (fTT)|) in two physically distinct scenarios: A the curvaton decays 
promptly into radiation and matter i.e. f-y(N) — 1, f m (N) — 1 and B curvaton decays first only into radiation 
component and the matter interaction begins when N — iV* i.e. f-y(N) = 1, f m (N) = 6(N — N*). In both scenarios 
we have systematically scanned the parameter space in order to identity physically acceptable parameter values. 

A. Continuous interactions, fj(N) = 1, f m (N) = 1 




-25 -20 -15 -10 -5 -25 -20 -15 -10 

log 10 (r 7 /ff ) Iog 10 (r r /i/ ) 

(c) (d) 



FIG. 1: Curvaton-radiation-matter system with different values of T 7 and F m , when (a) S7 CT o = 10 2 , (b) fi CT o = 10 5 , (c) 
flao = 10~ 7 and (d) Q^o = 10 -10 . Thick lines represent the initial system temperature in units of GeV. 

We have rigorously searched for initial values for which the system is physically motivated i.e. it passes the tests 



5 



mentioned in the previous section. This case has been previously studied in Q. In figures 1(a) |l(d) these tests can 



be seen as shades of black and white. A system that passes all of the above mentioned tests is labeled by white color 
whereas the opposite case is black. Shades of gray indicate that some but not all of the aforementioned physical 
requirements have been fulfilled. Besides these tests we have also included different contours of the initial system 
temperature into the figures. 

From the figures, we can read out that smaller values of initial curvaton density in general indicate higher initial 
temperatures, or higher reheat temperatures. The highest reheat temperatures correspond to higher values of F 7 , 
i.e. larger portion of the curvaton decays into radiation and therefore pushes the time of radiation-matter equality 
later leading to higher initial temperatures. Since the model applies only when the curvaton field is oscillating, this 
means that the temperature during reheating has been even higher. However as can be seen from figures, decreasing 
f2 CT an additional area becomes allowed corresponding to larger values of T m . This is a result of the fact that if T 7 
and T m are too small the curvaton does not decay fast enough, leading to possible issues during nucleosynthesis as 
mentioned above (e.g. for f2 CT o = 10~ 10 , T^/Ho = 10~ 20 and T m /Ho = 10~ 20 .) By increasing T m the curvaton field 
decays before nucleosynthesis and therefore leads to a physically sound system (e.g. for fl a o = ICP 10 , T 7 /Ho = 10~ 20 
and T m /H = 1(T 6 .) 

The figures indicate that if T 7im ~ H, a small initial curvaton density is required in order for the system to be 
physically viable. In addition, we see that in this case the reheat temperature is typically quite low. 

B. Time dependent interactions, / 7 (iV) = 1, fm(N) = 6(N - iV») 

When the curvaton-matter interaction is turned on at JV* the system is driven towards an equilibrium of (,„ — > 
—cf>(N^). This can be seen clearly from the equation of motion of i.e. eq. (fTT|) . because at times close to AT*, terms 
involving f' m (N) dominate the evolution and p a < 0. Thus the time dependent scenario resembles the situation we 
previously studied in [25|]. As can be seen from figures 2(a)||2(d) different initial values lead to allowed regions very 



similar to the previous case. Major alteration comes from the change of Hq — > H(N„) < Hq, which shifts different 
regions upward compared to the time independent scenario. Since the studied cases are initially radiation dominated, 
the Friedmann equation H 2 — p ~ p a + p-y = p a oe~ 3N + p 7 oe~ 4N ~ p 7 oe~ 4Ar allows us to estimate H(N st )/H ~ ICP 7 
when N* = 8. If the initial contribution of curvaton will be larger the system will become curvaton dominated earlier 
and therefore lead to a smaller value of H(N*). 

Another change can be seen in the initial system temperatures which rise steeply when T 7 increases. Still, higher 
values of T 7 lead to a larger contribution of radiation as in the time independent case. However, since now the 
matter interaction is turned on later, the curvaton might decay completely to radiation giving very low or no matter 
contribution at all. This in turn pushes initial temperatures up. Thus again we see that the initial curvaton density 
needs to be small. But a difference to the continuous-interaction case is that now the reheat temperature can be much 
higher. 

C. Isocurvature 

Isocurvature is defined as the difference between two curvature perturbations S mi = 3(C m — £ 7 ). Since the matter 
curvature perturbation Q m is here ultimately always driven to the value £ CT o [6J] the amount of isocurvature depends 
only on the final value of ^ 7 . The behavior of the radiation perturbation in the three- fluid model has been studied 
previously in Q . They found that if the curvaton begins to dominate before the decay epoch almost all of the radiation 
originates from the curvaton fluid and therefore gives £ 7 ~ Ca-o- 



We have plotted the amount of isocurvature in the three-fluid model in figures 3(a) and [3(b) Our results clearly 
agree with the reasoning above. If T m or T 7 is large enough, the curvaton fluid will decay before it begins to 
dominate the system and hence lead to a non-vanishing final isocurvature. For a smaller initial curvaton, it takes 
even longer for the curvaton fluid to dominate and therefore the system is adiabatic only for very small values of T m 
and T 7 . These in turn lead to high reheating temperatures, e.g. for initial values fl a o — ICP 10 , T^/Hq = ICP 26 and 

0.0006. The time 

dec 



F m /Ho — 10 34 a numerical evaluation gives a reheating temperature 10 12 GeV and S mj /( 



dependent interaction gives very similar results once the scaling T m /Ho — > T m /H* is taken into account. 



G 




FIG. 2: Curvaton-radiation-matter system where matter interaction is turned on at N — 8 with different values of F 7 and 
r m , when (a) Q, a o — 10~ 2 , (b) Q a o = 10~ 5 , (c) Slo-o = 10~ 7 and (d) fl a o = 10 -10 . Thick lines represent the initial system 
temperature in units of GeV. 



V. DISCUSSION AND CONCLUSIONS 



The curvaton model has gained a lot of attention in the recent years mainly because it can make the inflation 
potential look more natural . Since the first model where the curvaton decayed only into radiation [2l[ , a 
number of other possibilities have been explored including a curvaton web model the possibility of multiple 
curvaton fields [23| and different particle models such as axions [1^, [3(| just to name a few. 

In the present paper we have studied a three-fluid model in which the curvaton decays into both radiation and cold 
dark matter. This has been studied previously also in [f| [U [28j]. We have assumed that the initial system has no 
matter content and it is dominated by the radiation which originates from the decay of the inflaton field. Because all 
of the matter content of the universe comes from the curvaton field we are able to estimate the reheating temperature. 
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(a) (b) 
FIG. 3: Contours of the isocurvature/adiabatic ratio S m -y/Q = 3(Cm ~ £7)/?™ at the time of decoupling when (a) 

I dec I dec 

QcrO — 10~ 2 and (b) fla-o = 10~ 7 . In the darkest region the system is adiabatic whereas the white color corresponds to the 
largest amount of isocurvature. Larger initial contribution of curvaton clearly gives a larger adiabatic region. 



This can be additionally used to constraint the parameter space of the model. 

We have systematically scanned the parameter space and identified the regions where the model is physically 
acceptable, i.e. when evolution during and after nucleosynthesis is standard while requiring that the reheating 
temperature is not unreasonably high. We have identified these regions both when the decay rates are fixed and when 
the curvaton starts to decay later. These allowed regions are alike once the rescaling T/Hq — > T/H(N^) is taken into 
account and the real difference appears in the initial system temperature. 

We find that if the decay rates are comparable to the Hubble rate, a small initial curvaton density is required. 
Otherwise one needs to fine-tune the decay rates to be much smaller than H at the time of decay. In the continuous 
interaction case, requiring I\ ~ H leads to a low reheat temperature, but this can be avoided when the matter 
interaction is delayed. 

If the initial curvaton density is large, the final state is naturally adiabatic assuming that the system is otherwise 
physically acceptable. Note however, that this in turn requires fine-tuning in the decay rates. If I\ ~ H, we find that 
the final state generally contains a large isocurvature component. 

We have also we studied non-gaussianity in the framework of the three-fluid models. We find that in the region 
where the first-order perturbation theory can be applied, the three-fluid model gives no limits on the /jvl parameter. 
This is the result of a conserved curvature perturbation ( c , which carries the initial curvaton perturbation ( a into the 
matter perturbation £ m [fj. Our results differ from the previous results [1, H3] mainly because our f^L is evaluated 
at the time of last scattering and not at nucleosynthesis. This allows the non-gaussianity to be more easily compared 
to the observational Sachs- Wolfe effect and we do not have to use a radiation transfer function. In order to calculate 
the observational non-gaussianity, second-order perturbation theory needs to be applied [3l| . This is however beyond 
the scope of this article and will be the focus of a follow-up paper. 
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